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In the last years, stochastic Schrodinger equations have 
received considerable attention P, Q, M S IM 1^ H> IE 
IS ES in El El El El E3: two are the basic rea- 
sons for this increasing popularity. On the one hand, 
it has been shown (within the framework of dynami- 
cal reduction models) P, 0, H H ^3 that by adding 
commuting operators coupled to appropriate stochas- 
tic noises into the Schrodinger equation, the statevec- 
tor is driven into one of the common eigenmanifolds of 
such operators^; it is then possible to combine the quan- 
tum evolution and the process of wavepacket reduction 
into one single dynamical equation, providing a solution 
to the well known measurement problem of quantum 
mechanics. On the other hand, stochastic Schrodinger 
equations have proved to be a novel and useful tool 
in the study of open quantum systems, i.e. quantum 
systems interacting with the surrounding environment 
[illlllS,AJii,alElElElElE3 in this case, 
there are two advantages with respect to the standard 
formalism which resorts to the reduced density matrix of 
the open quantum system (when all the degrees of free- 
dom of the environment have been traced away): from a 
theoretical point of view, it is interesting to analyze how 
the statevector of such a system evolves in time, follow- 
ing the different possible realizations of the stochastic 
processes, while the reduced density matrix formalism 
gives only an ensemble description of the system. From 
a more "practical" point of view, in some cases it turns 
out to be simpler to unfold the evolution of a stochas- 
tic Schrodinger equation numerically on a computer and 
subsequently perform an average over sufficiently many 
different trajectories, rather than integrating numerically 
the equations for the corresponding statistical operator. 
Moreover, the statistical operator describing the ensem- 
ble of solutions of a stochastic Schrodinger equation is, 
by construction, a positive definite operator; this prop- 
erty is not always guaranteed by other approaches to the 



^ With an appropriate choice of the operators and noises, the re- 
duction of the statevector occurs only at the macroscopic level; 
the split between the micro and the macro is defined by the pa- 
rameters of the model. 



study of open quantum systems. 

In both stochastic Schrodinger equation takes 

the form 



/(L,z(0) 



\m), (1) 



where H is the free Hamiltonian of the system under 
study, while /(L, z(t)) is an appropriate function of the 
operators L = {Li, L2, . . . , L„} and of the stochastic 
noises z{t) = {21 (t), Z2{t), . . . , Zn{t)}. The statistical op- 
erator is defined as the ensemble mean, with respect to 
the noises z(i), of the projection operators \ip{t)){ip{t)\: 



pit) = {m)){m\)), 



(2) 



where (( 
equation 

d 



denotes the stochastic average. It obeys the 
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the form of the operator T[-] depending on the opera- 
tors L and the statistical properties of the noises z(t). 
In the case of statevector collapse models, T[-] describes 
the average effect of wavepacket reductions; in the case 
of open quantum systems, it embodies the effect of the 
surrounding environment on the system. 

In the white noise case, equations and |j2Jl take a 
simple form P H H 0, H S Hi ^-nd their most rele- 
vant properties have been extensively analyzed. How- 
ever, they do not encompass all physical interesting cases, 
since they describe only Markovian evolutions. 

Quite recently H, H, E El El El El El El , the 
formalism of stochastic Schrodinger equations has been 
generalized to the case of non-white complex Gaussian 
noises, which describe non-Markovian evolutions. The 
basic equation is: 



L • z(t) 



• / dssL{t,s) 



to 
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^ Note that, in general, it is not possible to write a closed equation 
for p{t). 



2 



where z(t) are complex gaussian noises with zero mean 
and correlation functions^: 

{{z*{t) z,{s))) = ay(i, s) {{zdt) z,{s))) = 0. (5) 

The dot • in equation Q denotes the product between 
the corresponding tensorial quantities, e.g. L • z(i) = 
T.^L^ z^{t). 

Equation I^J does not preserve the norm of the stat- 
evector and, as such, it does not have a direct physical 
meaning. It is possible to write an equation for the 
normalized physical vectors \(j){t)) = \ip{t)) /\\\tp{t))\\] this 
equation is rather involved since it does not have a closed 
form in \(f){t)), and is non linear. Anyway, for our pur- 
poses it suffices to work with Q), which gives the correct 
predictions once the probability distribution P[z(i)] as- 
sociated to the stochastic process of (O is replaced by: 



(6) 



Equation represents a breakthrough within the 
stochastic Schrodinger equations formalism, since it can 
describe a great variety of systems in many different phys- 
ical situations flit .13j . On the other hand, relations Q 
for the correlation functions set a limit for the possible 
choices of the noise, as they put on the same ground both 
the real and the imaginary parts of the complex noises"' 
z(i): this is evident once we re-write the correlation func- 
tions 10 in terms of x(t) = Re[z(t)] and y{t) = Im[z(t)]: 

{{x^{t)xj{s))) = {{yt{t)yj{s))) = ^Re[a^j{t,s)], 

{{x^{t)yj{s))) = -{{y,{t)xj{s))) = -^lm[a^j{t, s)]. 

It is then interesting to analyze what happens in the case 
of more general complex noises. There is a second reason 
why one would take into account more general stochastic 
Schrodinger equations, which is related to a characteristic 
feature of the formalism: in general, there are infinitely 
many different stochastic Schrodinger equations leading 
to the same equation for the statistical operator and, as 
such, reproducing the same statistical behavior. This is 
particularly simple to see in the white noise case, which 
we now briefly consider; for simplicity, we assume that 
the operators L are self-adjoint. In such a case, equation 
Q becomes^ jl]: 



dt 



\m) = 



L ■ z(t) - i 



(7) 



^ The complex functions a^j (t, s) are the coefficients of the matrix 

a(t,s) appearing in equation El . 
* As a matter of fact, the stochastic properties of z{t) are invariant 

under phase transformation, i.e. z(t) and z'{t) = e^'^z,(t) share 

the same statistical properties. 
® Note that all white noise equations are written as Stratonovich 

equations. 



where the correlation functions a{t, s) are now aij{t, s) = 
6ij S{t — s). The above equation, as previously mentioned, 
does not preserve the norm of the statevector; for com- 
pleteness, we write the corresponding norm-preserving 
equation for the physical vectors \(f>{t)): 



d 
di 



m)) = 



(L-(L)).(z(t) + (L)) 



:(L^-(L^)) 



(8) 



with (L) = {(l){t)\'L\(f){t)) . It is not difficult to prove that 
equation (jSJ or, equivalently, equation |7J plus the pre- 
scription ((nj, leads to the following equation for the sta- 
tistical operator^: 

|p(t) = -J[i7,p(t)]-i[L[L,p(i)]], (9) 

which of course is of the Lindblad type [2(Tj |. 

We now show that equation ^ can be derived from 
infinitely many different stochastic Schodinger equations. 
As a metter of fact, let us consider the noises w(t) = 
cx(i), where c = e"^ is a complex phase factor and x{t) 
is a real gaussian white noise with zero mean and cor- 
relation function {{xi{t) Xj{s))) — Sij6(t — s), and the 
following linear stochastic Schrodinger equation: 



-n" 



L • w(t) - cos^ if 



\m (10) 



A little calculation shows that, for any value of if, equa- 
tion (|10|l leads to Note that, in order to write equa- 
tion nni), we have resorted to complex white noises more 
general than those appearing in equation lO, since the 
correlation function {{wi{t) Wj{s))) = 6ij d{t — s) is no 
longer equal to zero (compare with equation lO). Note 
also that, in spite of the fact that they share the same 
statistical properties, equations (0) and ifTHI have a rad- 
ically different structure and lead to completely different 
time evolutions for the statevector^. 



® Th is is one of the basic equations in the theory of decoherence Il8l 
ITU; when the operators L commute among themselves, the non- 
Hamiltonian terms determine the damping of the off-diagonal 
elements (o|p(t)|/3), where \a), \f3) are two different eigenstastes 
of L. 

''' In refs. Il6l |21| the same issue is considered. In ref. |2 J . the 
authors prove that, in the white noise case, the two stochas- 
tic Schrodinger equations of the type 1101 with ip = (real 
noise) and 93 = 7r/2 (imaginary noise), lead to the same equa- 
tion for the statistical operator. In ref. Il6l the authors analyze 
two different unravelings for an open system interacting with a 
bath of harmonic oscillators, corresponding to two different sorts 
of measurements on the bath. The first unraveling leads to a 
non— Markovian stochastic Schrodinger equation of the type 
while the second unraveling leads to a non-Markovian stochastic 
Schrodinger equation with a real Gaussian noise (see equation 
1191 V By construction, the two equations describe (when the 
average is taken with respect to the noise) the evolution of the 
same quantum open system. 
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This "redundancy" in the number of stochastic 
Schrodinger equation leading to Q has an important 
consequence: one has the possibility to choose the most 
suitable equation according to the specific problem un- 
der study. For example, if one wants to reproduce the 
statevector reduction process, the best choice is equation 
(O , or equation (|10() with ip ^ (which corresponds to a 
purely real stochastic process): 



H + L-x(<) - 

h 



as already mentioned, the non-Hamiltonian terms of 
(fTT|l . together with the prescription ©, drive the stat- 
evector into one of the common eigenmanifolds of the 
operator L, if they commute among themselves. 

If on the other hand one wants to take advantage of 
the stochastic-Schrodinger-equation formalism to solve 
equation ||2J), then the mathematically simplest equation 
is I|1U|I with (y5 = 7r/2 (i.e. a purely imaginary noise): 



i L ■ x(t) 



\m)- (12) 



This equation is unitary, preserves the norm of the stat- 
evector and thus coincides with the physical equation for 
the normahzed vectors |0(<)} — |V'(i))/IIIV'(i))ll- The evo- 
lution it describes is different from that of equations l(TT|l 



and (UJ, since it does not lead to the localization of the 
statevector |21|. Anyway, as long as one is concerned 
only in time evolution of the density matrix, equation 
(112(1 is as good as (|llll and (in particular, confront 
the physical, norm-preserving Eqs. (|12l) and (jS))). 

In the light of the above discussion, it becomes interest- 
ing to generalize the stochastic Schrodinger equations for- 
malism to more general types of complex gaussian noises. 
To this purpose, let us consider the following stochastic 
Schrodinger equation: 



dt 



H + L-z(i) +0 



\m). (13) 



where, now, z(i) is a general complex gaussian stochastic 
process with zero mean and correlation functions: 



izt{t)z,{s))) = a,,{t,s) 
{{zi{t) Zj{s))) = bij{t,s); 



(14) 



aij{t,s) and bij{t,s) are generic functions, the only re- 
quirement on them being the consistency conditions: 
a^j{t,s) = a*^{s,t) and bij{t,s) = b^^{s,t). 

The operator O is uniquely defined 0, 0| by the re- 
quest that ifT^ preservers the average value of the square 
norm of \il){t)), this requirement being necessary in order 
to work out a meaningful physical interpretation: 



= mt)\ [L^ • ^*it) + ot + L . z{t) + o] m)))) 

= (((V'(01OlV'(0))) + |w01j[*rf^(Lt-a(i,s)+L.b(^,s)).^|VXt))^^ - 
mt)\0^ \m))) + |(VXi)l /* ds ■ (at(t, .) . L + bt(t, s) ■ Lt) \m) 

I 



0. (15) 



In the second line of equation H15|l , the functional deriva- 
tive acts on the ket \ip{t)), while in the third line it acts 
only on the bra ('ip{t)\. In going from the first to the sec- 
ond line, we have employed the Furutsu-Novikov formula 
[2^. Equation (|15|l is then satisfied if: 



O = - f ds {V ■ a(t, s)+L- h{t, s)) 

Jtn 



(5z(s)^ 



and equation (|13|) becomes: 



(16) 



We can easily see that equation (|17|l generalizes (@J and 
correctly reduces to it when b(t, s) = 0. One can write 
the equation for the evolution of the statistical operator 
and work out a perturbation expansion in more or less 
the same way it has been done for the b(i, s) = case 
[1^ , but this goes beyond the scope of the present 
article. 

It is instructive to analyze in greater detail the impor- 
tant case of self-adjoint operators L: 



dt 



m)) 



H + L- z{t) 



[ ds (Lt.a(i,s)+L-b(t,s)) 

J to 



(5z(s) 



m))- (17) 



dt 



\m) 



H + L-z(t) 



/ ds {a{t,s)+h(t,s)) 

Jtn 



(5z(s) 



\m). (18) 
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together with the hmit x(t) and y{t) 0; note that 
one cannot take such hmits for equation due to the 
particular form of the correlation functions 10 • In the 
first case (a purely imaginary noise), a{t, s) = —h{t, s), so 
that the third term at the right hand side disappears and 
(|18|l coincides with H12|) . This is an expected result, since 
a unitary evolution preserves the norm of the statevector 
and thus it does not need any corrective term. In the 
second case (a purely real noise), a(t, s) = h{t,s) and 
(fH^ becomes: 



dt 



\m) 



H 



L-z(t) 
II 

-2L- f dsa{t,s) 

Jtn 



Sz{s) 



This equation has been studied in detail in ref. [l5j, 
where it has been proved that — as for the white noise 
equation Ijllll — the non Hamiltonian terms drive the 
statevector into one of the common eigenmanifolds of the 
operators L, when they commute among themselves. 

To summarize, we have generalized the formalism of 
non-Markovian stochastic Schrodinger equations to any 
kind of coniplex Gaussian noises, recovering the results 
of 11,11 Hill 111 III III El as a special case. We have 
seen that such a formulation is more "flexible" as it covers 
the two particular but very important cases of purely real 
and purely imaginary noises, and in general it allows to 
take the most suitable form of the noise according to the 
specific problem under study. 



Electronic address: lbassi@ictp.trieste.it| 

G. C. Ghirardi, P. Pearle and A. Rimini, Phys. Rev. A 42, 

78 (1990). G.C. Ghirardi, R. Grassi and P. Pearle, Found. 

Phys. 20, 1271 (1990). G.C. Ghirardi, R. Grassi and A. 

Rimini, Phys. Rev. A 42, 1057 (1990). G.C. Ghirardi, R. 

Grassi and F. Benatti, Found. Phys. 25, 5 (1995). 

P. Pearle, Phys. Rev. Lett. 53, 1775 (1984). Phys. Rev. 

A 39, 2277 (1989). 

L. Diosi, Phys. Lett. A 132, 233 (1988). J. Phys. A 21, 
2885 (1988). Phys. Rev. A 42, 5086 (1990). 



[4] V.P. Belavkin, in Lecture Notes in Control and Infor- 
mation Science 121, A. Blaquiere ed., 245 (1988). V.P. 
Belavkin and P. Staszewski, Phys. Rev. A 45, 1347 
(1992). 

[5] A. Barchielli, Quantum Opt. 2, 423 (1990). Int. J. Theor. 
Phys. 32, 2221 (1993). 

[6] N. Gisin and J. Percival, Phys. Lett. A 167, 315 (1992). 
J. Phys. A 25, 5677 (1992). J. Phys. A 26, 2233 (1993); 
26, 2243 (1993). I. Percival, Quantum State Diffusion, 
Cambridge University Press, Cambridge (1998). 

[7] L.P. Hughston, Proc. R. Soc. London A 452, 953 (1996). 
S.L. Adier and L.P. Howritz: Journ. Math. Phys. 41, 
2485 (2000). S.L. Adler and T.A. Brun: Journ. Phys. 
A 34, 4797 (2001). S.L. Adler, Journ. Phys. A 35, 841 
(2002). D.C. Brody and L.P. Hughston, Proc. R. Soc. 
London A 458 (2002). 

[8] W.T. Strunz, Phys. Lett. A 224, 25 (1996). 

[9] L. Diosi, Quantum Semiclass. Opt. 8, 309 (1996). 
[10] L. Diosi and W.T. Strunz, Phys. Lett. A 2358, 569 

(1997) . 

[11] L. Diosi, N. Gisin and W.T. Strunz, Phys. Rev. A 58, 
1699 (1998). 

[12] W.T. Strunz, L. Diosi and N. Gisin, Phys. Rev. Lett. 82, 

1801 (1999). L. Diosi, N. Gisin, W.T. Strunz and T. Yu, 

Phys. Rev. Lett. 83, 4909 (1999). 
[13] T. Yu, L. Diosi, N. Gisin and W.T. Strunz, Phys. Rev. 

A 60, 91 (1999). Phys. Lett. A 265, 331 (2000). 
[14] A. A. Budini, Phys. Rev. A 63, 012106 (2000). 
[15] A. Bassi and G.C. Ghirardi, Phys. Rev. A 65, 042114 

(2002). 

[16] J. Gambetta and H.M. Wiseman, Phys. Rev. A 66, 
012108 (2002). J. Gambetta and H.M. Wiseman, to ap- 
pear in Phys. Rev. A (preprint quant-ph/0208169). 

[17] M.B. Plenio and P.L. Knight, Rev. Mod. Phys. 70, 101 

(1998) . 

[18] E. Joos and H.D. Zeh, Zeit ffir Phys. B 59, 233 (1985). 

[19] D. Giulini, E. Joos, C. Kiefer, J. Kupsch, I.-O. Sta- 
matescu and H.D. Zeh, Decoherence and the Appear- 
ance of a Classical World in Quantum Theory, Springer 
(1996). 

[20] G. Lindblad, Commun. Math. Phys. 48, 119 (1976). 
[21] G.C. Ghirardi and R. Grassi, in: Nuovi problemi della 

logica e della filosofia della scienza, CLUEB, Bologna 

(1991). 

[22] A. Novikov, Zh. Eksp. Teor. Fiz. 47, 1915 (1964) [Sov. 
Phys. JETP 20, 1290 (1965)]. K. Sobczyk, Stochastic 
Wave Propagation, Elsevier, Amsterdam 1985. 



